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OO I ABSTRACT 

j> ■ We present 350 //m observations of 36 ultraluminous infrared galaxies 

^ ■ (ULIRGs) at intermediate redshifts (0.089 < z < 0.926) using the Submillimeter 

■ High Angular Resolution Camera II (SHARC-II) on the Caltech Submillimeter 
^ ■ Observatory (CSO). In total, 28 sources are detected at S/N > 3, providing the 

■ first fiux measurements longward of 100 /xm for a statistically significant sam- 
^ ■ pie of ULIRCs in the redshift range of0.1<z<1.0. Combining our 350 /xm 

P^. fiux measurements with the existing IRAS 60 and 100 /xm data, we fit a single- 

Q . temperature model to the spectral energy distribution (SED), and thereby esti- 

. mate dust temperatures and far-IR luminosities. Assuming an emissivity index 

. of /? = 1.5, we find a median dust temperature and far-IR luminosity of = 

> : 42.8±7.1K and log(LFiR/L0) = 12.2±0.5, respectively. The far-IR/radio cor- 

^ . relation observed in local star-forming galaxies is found to hold for ULIRCs in 

the redshift range 0.1 < z < 0.5, suggesting that the dust in these sources is 
predominantly heated by starbursts. We compare the far-IR luminosities and 
dust temperatures derived for dusty galaxy samples at low and high redshifts 
with our sample of ULIRCs at intermediate redshift. A general Lfi^-T^ relation 
is observed, albeit with significant scatter, due to differing selection effects and 
variations in dust mass and grain properties. The relatively high dust tempera- 
tures observed for our sample compared to that of high-z submillimeter-selected 
starbursts with similar far-IR luminosities suggest that the dominant star forma- 
tion in ULIRCs at moderate redshifts takes place on smaller spatial scales than 
at higher redshifts. 



Subject headings: galaxies: starbursts — infrared: galaxies — infrared: ISM: 
continuum 
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1. Introduction 

Dust grains in the interstellar medium (ISM) of galaxies are heated by various energy 
sources in vastly different physical environments: the diffuse interstellar radiation field, active 
star formation and active galactic nuclei (AGN). The bulk of the emission from heated 
dust grains occurs at far-IR and submillimeter (submm) wavelengths, and can account for 
a significant, if not dominant, fraction of the total energy output from dust-enshrouded 
galaxies where extinction and re-radiation by dust lead to high IR-to-optical flux ratios. 

In the latter category we find the population of ULIRGs (log (Lj^/L©) > 12), which 
were first discovered in large numbers in the local Universe by IRAS (Soifer et al. 1986). 
These galaxies harbor vast amounts of dust, responsible for obscuring most of the intense 
star formation and AGN activity that are taking place. Observations in the optical/near-IR 
as well as of rotational CO line emission revealed that ULIRGs tend to be associated with 
interacting or merging gas-rich disks, with the more luminous being the more disturbed and 
morphologically complex systems (Sanders & Mirabel 1996), thus suggesting that ULIRGs 
are driven by mergers in which gas is funneled towards the central regions, fueling an intense 
starburst and/or a black hole. While local {z < 0.1) ULIRGs are believed to re-enact some 
of the processes relevant for galaxy formation at high redshifts, it is essential that studies of 
ULIRGs over a wide redshift range are carried out, in order to get a complete picture of how 
the most massive galaxies formed as a function of redshift. At high redshifts (1 < 2; < 3), 
the population of optically faint radio sources - which includes (sub)mm-selected galaxies 
(SMGs) uncovered by SCUBA and mambo (Small et al. 1997; Hughes et al. 1998; Bertoldi et 
al. 2000; Carilli et al. 2001) as well as the population of optically faint 24 //m sources selected 
with Spitzer (Papovich et al. 2004) - is thought to represent the high-z analogues of local 
ULIRGs. 

In contrast to the populations of local and distant ULIRGs, the population at intermedi- 
ate redshifts (0.1 < z < 1.0) remains virtually unexplored at far-IR/submm wavelengths. This 
intermediate redshift interval is of particular interest as it bridges the gap between the distant 
and local ULIRG populations, and broadly marks the transition period where ULIRGs cede 
their position to the less extreme luminous infrared galaxies (LIRGs - log (Li^/L©) > 11) as 
the dominating contributor to the co-moving IR luminosity density (Le Floc'h et al. 2005). 
In addition, the cosmic infrared background (CIRB - Hauser et al. 1998; Schlegel et al. 
1998), which peaks at ^150 /im and constitutes about half of the total energy density in 
the Universe, is dominated by emission from (U) LIRGs in the above intermediate redshifts 
range (Lagache et al. 2005). 

In this paper we present 350 fim observations of a statistically significant sample of 36 
ULIRGs in the intermediate redshift range 0.089 < z < 0.926. Throughout this paper we 
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adopt a flat cosmology with Hq = 70kms ^ Mpc ^, = 0.30 and Q\ = 0.70. 

2. The Sample 

Our sample is a subset of the FSC-FIRST (FF) catalog (Stanford et al. 2000), which was 
compiled by cross-correlation between the IRAS Faint Source Catalog (FSC - Moshir et al. 
1992) and the Faint Images of the Radio Sky at Twenty cm (FIRST - Becker et al. 1995), and 
consists of 108 ULIRGs for which spectroscopic redshift information and X-band imaging 
were successfully obtained. Situated in the redshift range of 0.1 < z < 1, the FF sources 
constitute a well deflned and statistically signiflcant sample of ULIRGs at intermediate 
redshifts. Only sources with secure IRAS detections at both 60 and 100 /xm were considered 
for 350 /xm observations, with no flux cutoflF employed at either waveband. The fact that 
virtually no other far-IR/submm continuum data points exist longward of 100 fim for the FF 
sources, makes our 350 /xm flux measurements extremely valuable for attaining meaningful 
estimates of the SEDs. Furthermore, the availability of the spectroscopic redshift information 
allows accurate estimates of dust temperatures and far-IR luminosities. Finally, we note 
that targeted observations at 350 /xm were made possible due to the positional accuracy 
provided by the existing radio observations. This is highly desirable for far-IR/submm 
observations given the limited angular resolution currently available to telescopes working 
in this wavelength range. In total, we observed 36 sources in the FF catalog at 350 /xm, 
which are listed in Table [U along with their redshifts and flux densities at 60 /xm, 100 /xm 
and 1.4 GHz. The redshift distribution of our sample, having a sample median of 0.366, is 
shown in Fig. [H 

There are a few selection biases possibly present in the FF catalog. First, high-redshift 
galaxies are potentially underrepresented due to the adverse X-correction at 60 /xm and 
1.4 GHz. Second, since the FSC is a flux-limited sample (the 60 /xm cutoff flux density is 
^ 0.2 Jy - Moshir et al. 1992), the selection function is sensitive to dust temperature and 
biased against galaxies with low dust temperatures. Third, galaxies with enhanced AGN 
activity are probably overrepresented in the sample, since such sources are expected to be 
more luminous (Genzel et al. 2000). 

3. Observations 

The observations were made during a series of SHARC-II (Dowell et al. 2003) observing 
runs from January 2003 to September 2004 at the CSO under excellent weather conditions 
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("^225 GHz ^ 0.06). Integration time varied from source to source, depending on the brightness 
of each source at 350 /xm, atmospheric opacity and sky variabihty. On average, 1.5 hours 
of integration time was spent on each source. AU measurements were made by scanning 
the bolometer array in a Lissajous pattern centered on each source. With the exception of 
FF 0738+2856, FF 0758+2851, and FF 0804+3919, scans were taken without chopping the 
secondary mirror. Pointing and cahbration scans were taken on an hourly basis on strong, 
point-like submillimeter sources. SHARC-II raw data - taken on both calibration and science 
sources - were reduced by using the software package CRUSH (Kovacs 2006), with the same 
''crush" option "-deep". Every single calibration scan was reduced individually, to obtain 
a 350 /xm map and flux density. The instrument flux was compared to the true flu:x§, and 
a calibration factor, deflned as X^^i = S'true/^'mstrument, was calculated. We also applied the 
built-in CRUSH utility "calibrate" to the 350 /xm calibrator map, to have a measure of the 
actual beam size (FWHM). All X factors and FWHM values thereby estimated were complied 
for each evening, from which calibration scans that yielded 0.5 < X < 2 and FWHM < 20'^ 
were deemed reliable, and thus selected to form a subset. Finally, the sample median of 
the Xcai factor subset is taken as the Xcai factor appropriate for use to reduce science scans 
taken on the same evening. The calibration accuracy, as measured by the standard deviation 
of the Xcai factor subset, was found to be generally within 20%. Consistent calibrated flux 
measurements were also found between the chopped and unchopped observations. 

In total, 28 out of the 36 sources in our sample were detected with signal-to-noise ratios 
S/N > 3. This high detection rate (^-^78%) demonstrates the effectiveness of SHARC-H 
for studying faint 350 /xm sources at the 50 — 100 mJy level. Other than the IRAS data, 
the present results at 350 fim are the flrst reported detections of FF sources at submm 
wavelengths longward of 100 /xm, except for FF 1532+3242 and FF 1614+3234 which were 
previously observed at 90 and 180 jam by Verma et al. (2002), and at 850 /xm by Farrah et al. 
(2002) respectively. Table [1] lists the 350 /xm flux measurements of the FF sources selected 
for SHARC-II observations. The final 350 /xm signal maps smoothed with a 9^^ FWHM beam 
are shown in Figs. [2] and [31 

4. Spectral Fits and Derived Quantities 

We adopt the single-temperature, optically-thin far-IR/submm SED model (Hildebrand 
1983) 

S. = f^B.(Td)Q., (1) 



^Available from ^http: / / www.submm.caltech.edu/sharc/ 
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where is the dust temperature, = Qo (— )^ is the absorption coefficient and f] is the 
emissivity index. A flux error at a uniform 20% level is assumed at each data point E, 
hence each flux measurement is weighted equally in the nonlinear least squares flt. It is not 
possible to flt both and /? as free parameters due to the paucity of available data points. 
Additionally, the negative T^-/3 degeneracy in the SED fltting procedures (Blain et al. 2003) 
often renders simultaneous estimates of and (3 problematic. To achieve redundancy, we 
flt three SED curves for each source in our sample with /? flxed to 1.0, 1.5 and 2.0. The 
fltted SEDs for each source are overplotted (Fig. Hj) whenever convergence can be achieved 
by the fltting procedure. (3 values favored by the observed photometric data points diflFer 
from source to source, suggesting real diflFerences in emissivity index among the galaxies. 
The negative T^-(3 degeneracy is also clearly manifested in Fig. [H which shows that higher 
estimates of T^ are always associated with lower values assumed for /?. We assume /? = 1.5 - 
a value most commonly assumed and observed - to estimate eflPective dust temperature and 
derive physical quantities in later analysis and discussion (Table [1]). 

The best-fltting dust temperatures in the SHARC-II detected FF sample, with /? assumed 
to be 1.5, are found to range within a wide band of 35.1 < T^ < 62.6 K, with a sample me- 
dian and standard deviation T^ = 42.8 ± 7.1 K. A histogram of T^ is shown in Fig. [5l 
All SHARC-II detected FF galaxies have T^ values in the range of 35 to 50 K, except for 
four sources - FF 0050-0039 (z=0.727), FF 1042+3231 (z=0.633), FF 1106+3201 (z=0.900), 
FF 1532+3242 (z=0.926) - that have Td > 50 K. The apparent absence of low-Td-high-z 
galaxies within the sample could at least in part be caused by the selection eflFects inherent 
to the FF catalog, while the lack of high-T^-low-z (effectively high-Td-low-LpiR) galaxies in 
the sample most likely reflects real properties of the ULIRG population. 

The global far-IR luminosities (Lfj^) and dust masses (M^) are calculated using 

/^lOOO /im 

= AnDl / S^du (2) 

where Dl is the luminosity distance, and k = is the dust mass absorption coeflicient, 

varying as oc u^. We adopt dust parameters given by Hildebrand (1983): Qu^j^m = 7.5x 10~^, 
p = 3.0g/cm^ and a = O.l/xm, which implies = 1.875 (kg/m^)~-^. While there 

are signiflcant uncertainties associated with due to the lack of accurate knowledge of 
interstellar dust properties, the inferred far-IR luminosities are tightly constrained as long 



^Although the IRAS 100 /im measurement are probably less reliable than the 60 and 350 /im data points. 
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as a satisfactory SED fit is achieved. The Lpm and values derived for our sample lie in 
the ranges 10^^-^ < Lpu, < 10^^-^ and lO^-^ < Md < 10 ^-^Mq, respectively, with sample 
medians of Lpi^ = 10 ^^-^^^-^L© and = 10 ^-^^^-^ M©. Thus, our sample spans nearly two 
orders of magnitude in far-IR luminosity, and more than an order of magnitude in dust mass. 

If the far-IR luminosity originates from warm dust heated by massive, short-lived stars, 
which we argue in §5.11 is the case for majority of our sources, it provides a measure of the 
instantaneous star formation rate (SFR) in these galaxies. Adopting the calibration factor 
given by Kennicutt (1998), we have 

= 1 7 X ^^^^ fX U) 

where L^m is integrated over the wavelength range of 40 — 1000 /xm, / characterizes the 
percentage energy contribution to far-IR luminosity by massive young stars (/ < 1), and 
X accounts for integrated fiux shortward of 40 /xm {X > 1). Using the estimated value 
of X ^ 1.4 (from Eqs. [3 and [8]), we find the median star formation rate of our sample 
is SFR = lO^'^^^"^ M0 yr~^. The implied star formation rates are thus 3 — 4 orders of 
magnitude higher than what is typically observed in in normal quiescent galaxies (Kennicutt 
1998), but comparable to that of the local ULIRG Arp220 (SFR ^ 320M yr"^ - Dopita 
et al. 2005). Such intense starbursts can be induced by strong tidal interactions, a scenario 
consistent with high-resolution radio and near-IR observations of the FF catalog that reveal 
morphological features typical of strong interacting/merging systems for the majority of the 
sources (Stanford et al. 2000). 



5. Discussion 

5.1. Far-IR/Radio Correlation 

A tight far-IR/radio correlation is known to hold for a wide range of star- forming galaxies 
in the local Universe (Condon 1992). This is explained by the fact that radio emission at 
1.4 GHz is thought to arise mainly from nonthermal synchrotron emission of cosmic electrons 
that are released by supernovae explosions of massive young stars, thus radio emission is 
expected to be proportional to the present star formation rate, as is far-IR emission. Whether 
and how this far-IR/radio correlation might evolve with galaxy characteristics and redshift 
remains an open question at the present, partly due to the substantial observational efforts 
involved to compile the relevant information (i.e., multiband far-IR/submm fluxes, 1.4 GHz 
radio flux, and precise spectroscopic redshifts) for a large number of galaxies at intermediate 
and high redshifts. Our sample of SHARC-II detected FF sources, having the complete 
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information mentioned above, constitutes a statistically significant sample that is suitable for 
a meaningful investigation of the far-IR/radio correlation at moderate redshifts 0.1 < z < 1.0. 

The far-IR/radio correlation for our sample is plotted in Fig. [6^. Note that the corre- 
lation is investigated over the full far-IR/submm wavelength range of 40 — 1000 /xm. The 
far-IR/radio correlation is often expressed in terms of the logarithmic far-IR/radio flux ratio 
q (Helou et al. 1985), deflned as 

where Li,4ghz is the 1.4 GHz (rest-frame) radio luminosity. Accounting for both the K- 
correction, through (1 + z)^, and the bandwidth compression by a factor of (1 + z)~^ ^ we 
have 

'°8< WlP> = 1^ 0« + 2 log (^) + („ - 1) log (1 + z) + log(5HSii) , (6) 

where a, deflned by S^^ oc is the radio spectral index. Realistic constraints on a 

are not currently available for the FF sample due to the lack of multi-wavelength radio 
measurements. We therefore assume a constant value oi a ^ 0.7 in the spectral domain 
near 1.4 GHz - a value that is appropriate when synchrotron emission is dominant (Condon 
1992). Following Eqs. ([5j) and ([6j), and using rest-frame L4o_iooo/im, we have derived q for 
each source in our sample of SHARC-II detected ULIRGs (Tabled]). A broad range in q 
values are found (1.44 < q < 2.76), with a sample median oi q = 2.4 ± 0.4. The distribution 
of q values for our sample is shown in Fig. (6)3. Of the entire sample, 21 sources (75%) have q 
values within ±lcr of the median value, while the seven remaining outliers all have q < 2.0, 
suggesting that they might be radio-loud AGNs. We note that four of the seven outliers 
are at z > 0.5 and flve have Lpi^ > 10^^-^ L©, which is most likely caused by the selection 
biases in our sample that favor radio-loud, luminous objects with enhanced AGN activity 
at high redshifts (^. The suspected AGN fraction of ~25% is also roughtly consistent 
with the ^^10% found for the entire FF sample, and with the trend of an increasing AGN 
fraction with FIR-luminosity as noted by Stanford et al. (2000). We therefore exclude the 
seven sources in the investigation of the far-IR/radio correlation by setting a lower bound 
of (g) — = 2.0. To check this approach we have searched the literature for independent 
AGN/starburst classiflcation for our sources. In total, classiflcations are found for six sources 
(Stanford et al. 2000; Zakamska et al. 2004; Best et al. 2005; Best private communications), 
of which three are classifled as AGNs (FF 0050-0039, FF 0804+3919 and FF 1532+3242) and 
three as starbursts (FF 0030-0027, FF 0907+3931 and FF 2221-0042). The three known 
AGNs have q < 2.0, while the three known starbursts are within ±lcr of our median q value, 
in agreement with our lower cut in q. However, we wish to note that a systematic change in 
the radio slope would have a direct impact on the estimated q values. For instance, Ivison et 
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al. (private communications) have suggested that SMGs might have shaUower synchrotron 
slopes than "non-SMG" radio sources. In this scenario, q values estimated by assuming 
a ^ 0.7 would be underestimated, with the underestimates in q being more severe at higher 
redshifts. As such, the very low q values {q < 2.0) characteristic of the high-z {z > 0.5) 
sources in our sample could also be caused by actual radio slopes that are flatter than the 
assumed value of 0.7 in these sources. 

The remaining 21 sources with q > 2.0 show a very tight far-IR/radio correlation (Fig. 
[Gb), with a sample median of q = 2.6 ±0.2. We wish to compare the far-IR/radio correlation 
observed at intermediate redshifts to that in the local Universe, for which q values are found 
to be bound within a narrow range 2.3 ± 0.2, with the far-IR luminosity deflned to cover the 
wavelength range of 40 — 120 /xm (Condon et al. 1991). To account for the difference in the 
wavelength integration range, we adopt estimates given by Soifer & Neugebauer (1991) 

L8-40/im ^ 40% X L40_i000/im (7) 
Ls-lOOO/im ^ (1 + 40%) X Ls-m^m, (8) 

from which we derive ^ 

L40-1000/im ^ — X L4o_i20/im- (9) 

The local far-IR/radio correlation as given by Condon et al. (1991) can thus be recast into q 
values with a median (q) = 2.3 + log {■^) = 2.5 and scatter = 0.2, in excellent agreement 
with the q values calculated for our sample of intermediate redshift ULIRGs. Furthermore, 
a flt to log (Lfir) and log (Li.4ghz) inferred for our sample yields a slope of a = 1.00 ± 0.14, 
which is also consistent with that of the local correlation (a = 1.11 ± 0.02 - Condon et al. 
1991). 

Thus, we conclude that our derived far-IR/radio correlation in the redshift range 0.1 < 
z < 0.5 matches very well with that of the present day, suggesting that the correlation 
does not evolve signiflcantly from z ^ to z ^ 0.5. This supports the notion that dust 
heating in the intermediate-redshift ULIRGs originates predominantly from star formation, 
with modest contribution to the far-IR luminosity from old stellar populations and AGNs. 
Unfortunately, we cannot confldently extend this conclusion to z ^ 1, since four out of the 
five sources in our sample at z > 0.5 were discarded on suspicion of being radio-loud AGNs, 
albeit the possibility of flatter radio slopes in the high-z sources. Second, we flnd no evidence 
of an increase in the scatter of the correlation with redshift. Our study validates the approach 
by Appleton et al. (2004), in which Spitzer 70 /xm data are X-corrected in a model-dependent 
manner to derive the full far-IR luminosities for sources out to z ^ 1. Such conflrmation is 
important, allowing further speculation with a single flux density measurement. 
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5.2. Luminosity- Temperature Relation 

The far-IR luminosities and dust temperatures inferred for our sample are plotted 
against each other in Fig. [71 For comparison, we have also plotted Lpi^ and T^ values 
compiled by Blain et al. (2004) for samples of luminous infrared galaxies at low and high 
redshifts. In the case of the low-z sample, 60, 100 and 850 /xm data are available and utilized 
in the SED fitting (Dunne et al. 2000). By contrast, Lpm and T^ values of the high-redshift 
sample are derived by assuming the far-IR/radio correlation (Chapman et al. 2003; Chapman 
et al. 2005), which has been subsequently confirmed by Kovacs et al. (2006). 

Taken at face value, each sample exhibits an overall trend between Lpi^ and T^, where 
sources with large far-IR luminosities having, in general, higher dust temperatures. This 
trend can be understood in terms of the approximate relation 

L™oc«:oMdT^+^ (10) 

where /^o is the dust mass absorption coefficient normalized at some reference frequency 
However, there is significant scatter in the Lpi^-Td relation, due to variations in the dust 
emissivity, as well as in the amount of dust in each galaxy. 

Although significant scatter in T^ for a given Lpm exists for galaxies within each sam- 
ple. Fig. [3 reveals systematic shifts between the three samples shown. In terms of dust 
temperature, the low-z sample spans a fairly narrow range (T^ ^ 25 — 45 K) with a median 
temperature of T^ = 35.6 ± 4.9 K (Dunne et al. 2000). The latter is consistent with the me- 
dian temperatures found for our sample (T^ = 42.8 ±7.1 K) and the SMGs (T^ = 36.0 ±7.0 K 
- Chapman et al. 2005), although the latter two span a larger range of temperatures. In 
particular, we notice that all sources in our sample have T^ > 30 K, whereas ^30% of the 
SMGs have dust temperatures below 30 K. Furthermore, there are four sources in our sample 
for which we can confidently say they have dust temperatures above 50 K, while all SMGs 
are consistent with having T^ < 50 K given their large error bars. Thus, it seems that our 
sample overall consists of sources with somewhat warmer dust than both the local and high-z 
sample. Assuming the dust emissivity properties in starburst galaxies do not change sub- 
stantially as a function of redshift, and by Eq. (fTOl) . SMGs thus appear to contain cooler but 
considerably larger amounts of dust than (U)LIRGs at low and moderate redshifts (Fig. [Tj). 



' LpiR = ^^^oMd(kTd/h^o)^+^r(4 + /3)C(4 + /3), where r(z) = t^-^e-^t and 
^(^) ~ r(s) lo^ Gamma and Riemann ( functions respectively. This formula is equivalent 

to that given by De Breuck et al. (2003). With the normal range of 1.0 < < 2.0, C(4 + (3) is nearly flat 
while r(4 + /3) varies rapidly. Therefore the approximation as stated in Eq. (p!Q|) is valid when variations in 
/3 are negligible. 
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This is consistent with observations of CO rotational line emission in bright SMGs, which 
have revealed that they are very massive systems, having at least four time more molecular 
gas than that in local ULIRGs (Neri et al. 2003; Greve et al. 2005). We argue the observed 
difference can be explained as the combined effects of selection effects and intrinsic evolution 
in the ULIRG population as a function of redshift. 

First, the very high far-IR luminosities characteristic of the high- 2; SMG sample can 
be caused by the flux cut-offs of present-day submm and radio surveys, which limit the 
detections of high-z sources with Lpm < 10^^ L©. Second, by virtue of their selection at 
submm wavelengths, the SMGs are naturally biased in favor of cooler dust temperatures. 
This has been confirmed by recent Spitzer imaging, which shows that the vast majority 
80 — 90%) of the Spitzer detected SMGs have near-IR/mid-IR properties typical of cold, 
starburst-like ULIRGs (Egami et al. 2004). However, it has been suggested that a significant 
population of optically faint radio galaxies (OFRGs), which so far have remained undetected 
by current 850 fim surveys due to their warmer dust (T^ > 50 K), might exist and constitute 
an extension of SMGs to warmer temperatures (Chapman et al. 2004b). OFRGs are luminous 
in the radio and the 24 fim wavebands, and it seems our intermediate redshift sample, being 
select by IRAS^ might be more representative of OFRGs than SMGs. Finally, the existence 
of low-to-intermediate redshift (U)LIRGs with cool dust temperatures cannot be ruled out, 
since such sources are unlikely to have made it into the 60 and 100 /im /i?^S'-selected samples 
of local (U)LIRGs. The selection biases are more severe for the intermediate-redshift sample 
than for the local sample, as predicted by the shape of the iC-correction at these wavebands. 
Unbiased wide surveys with SCUBA-II and Herschel at 200 — 850 /im are expected to uncover 
cool sources at z ^ — 1 if they exist, thus allowing for a more direct comparison with the 
high- 2; population. 

While at least part of the observed differences in dust masses and temperatures between 
the various samples can be attributed to the selection effects specific to each sample, real 
cosmological evolution of the ULIRG population with redshift appears to play a significant 
role. This is strongly evidenced by the lack of very luminous (Lpm > 10-^^ Lq) sources in 
the \ow-z sample, which reflects the well-known decrease in the abundance of such extreme 
systems by about three orders of magnitude from z ^ 2 — 3 to the present day. At moderate 
redshifts, such systems are more common as suggested by the overlap in far-IR luminosity 
between our sample and the high-z SMGs. 



5.3. Dust Temperature and Physical Scale of Star Formation 



Finally, we wish to argue that the higher dust temperatures observed in our sample of 
ULIRGs at moderate redshifts, as compared to those of the high-z SMGs, may be linked 
to variations in the characteristic physical scale on which star formation dominates at these 
two different epochs. Eq. (fTOl) suggests that the observed T^, characterizing the large-scale 
SED of a galaxy is determined by the total dust heating (and by implication the global star 
formation rate) per unit dust mass within the galaxy The dust temperature is therefore 
linked to the global star formation efficiency subject to uncertainties in dust emissivity 
additional dust heating, and gas-to-dust ratio. At the same time, the star formation efficiency 
is known to be primarily regulated by gas density in star-forming regions (Schmidt law) across 
enormous ranges spanning 5 — 6 orders of magnitude in gas densities and star formation 
rates (Kennicutt 1998). Calculations by Yang (2006) have shown that a general relationship 
can be expected to exist between the observed T^ and the spatial extent of star formation, 
suggesting higher observed T^ are generally indicative of star formation that occurs over more 
concentrated regions. This proposed link is in good agreement with observational evidence; 
high-resolution radio observations of (cooler) high-redshift SMGs show spatially extended 
starbursts on scales of ^lOkpc (Chapman et al. 2004a), while the few local ULIRGs with far- 
IR luminosities comparable to that of SMGs in general have higher dust temperatures with 
the dust and molecular gas being concentrated in compact regions with scales on the order 
of 0.1^1 kpc (Young & Scoville 1991; Sanders & Mirabel 1996). Today's mm interferometers 
such as the IRAM Plateau de Bure Interferometer (PdBI) have just about the sufficient 
sensitivity as well as spatial resolution (FWHM:^ 0.6" at 1.3mm for PdBIin its extended B 
and A configuration) to allow us to detect and marginally resolve the dust emission in sources 
such as the ones in our sample. However, such observations are extremely challenging and 
time consuming, and will at best provide us with only a rough gauge, or even perhaps only 
an upper limit, of the physical extent of the dust. Of course, once ALMA comes on-line 
the situation will improve dramatically, resulting in detailed maps of the dust emission that 
will not only pinpoint the hot dust, which presumably originates from compact regions of 
intense star formation and/or AGN- activity, but also outline the morphology and scale of 
any underlying, cooler dust reservoir that may exist. 



6. Conclusion 

This paper reports 350 /xm observations of 36 intermediate-redshift ULIRGs (0.089 < 
z < 0.926), out of which 28 are detected. The newly acquired 350 /xm data are the first 
reported flux measurements longward of 100 /xm for these galaxies, and lead to mean- 
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ingful estimates of their dust temperatures (T^ = 42.8 ± 7.1 K) and far-IR luminosities 
(LpiK = lO^^*^^^*^ L0). These intermediate-redshift ULIRGs appear to host large amount 
of warm dust (M^ = 10 ^-^^^-^ M©) and are experiencing intense starburst activity (SFR = 
2q2.5±o.5 yr~^). The far-IR/radio correlation observed in star-forming galaxies in the local 
Universe remains valid for ULIRGs in the more distant Universe over the redshift range of 
0.1 < z < 0.5, where dust heating seems predominantly powered by ongoing star formation, 
with modest contribution to the far-IR luminosity from older stars and AGNs. Far-IR lu- 
minosities and dust temperatures derived for dusty galaxies over a wide range of redshifts 
reveal a general positive relation between theses two quantities. However, there is significant 
scatter in the Lfir-T^ relation, due to differing selection biases as well as variations in dust 
mass and grain properties. We argue that the somewhat higher dust temperatures of our 
sample compared to the high-z SMGs are possibly tied to the more compact spatial scales 
of the ongoing star formation. 

We thank Attila Kovacs for help with data reduction and taking data on our behalf; 
Andrew Blain for useful discussions and providing part of the Lpm-Td data; and Philip Best 
for providing us with AGN/starburst classifications for some of our sources. We are deeply 
grateful to the referee for insightful and thorough comments. The CSO is supported by the 
NSF fund under contract AST 02-29008. 
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Table 1. FF sources selected for SHARC-II observations (36). 



Source Name RA Dec z Seo/xm"" Sioo/xm^ Ssso/xm Si.4ghz log LpiR log q 

J2000 J2000 Jy Jy mJy mJy K L© 



FF 0030-0027 


00 


30 


09, 


.099 


-00 


27 


44. 


,40 


0, 


.242 


0, 


,59 


1, 


,11 


134 ± 16.8 


2, 


,52 


40, 


,4 


± 


2, 


.1 


12 


.25 


8.28 


2, 


.64 


FF 0050-0039 


00 


50 


09, 


.806 


-00 


39 


00. 


,96 


0, 


.727 


0, 


,22 


0. 


,49 


43.0 ± 10.4 


4, 


,32 


57. 


,6 


± 


2, 


.8 


12, 


.84 


8.15 


1, 


.91 


FF 0123+0114 


01 


23 


06, 


.973 


+01 


14 


10. 


,03 


0, 


.089 


0, 


,23 


0, 


,55 


( ±9.75) 


1, 


.78 






















FF 0240-0042 


02 


40 


08, 


.576 


-00 


42 


03. 


,56 


0, 


.410 


0, 


,20 


0, 


,69 


66.5 + 10.4 


0. 


,84 


42, 


,9 


± 


1, 


.8 


12, 


.40 


8.30 


2, 


.75 


FF 0245+0123 


02 


45 


55, 


.355 


+01 


23 


28. 


,40 


0, 


.798 


0, 


,17 


0, 


,32 


(32.4 + 25.6) 


2. 


,03 






















FF 0312+0058 


03 


12 


38, 


.445 


+00 


58 


33. 


,86 


0, 


.130 


0, 


,26 


0, 


,61 


(19.2 + 12.1) 


3, 


,84 






















FF 0317-0129 


03 


17 


43, 


.635 


-01 


29 


07. 


.33 


0, 


.265 


0, 


.21 


0, 


.34 


106 ± 9.84 


1, 


,67 


36, 


,2 


± 


2, 


.2 


11 


.91 


8.18 


2, 


.39 


FF 0738+2856 


07 


38 


29, 


.856 


+28 


56 


38. 


,74 


0, 


.334 


0, 


,25 


0, 


,64 


68.7 + 8.41 


1, 


.38 


41, 


.9 


± 


1, 


.9 


12, 


.24 


8.20 


2, 


.59 


FF 0748+3343 


07 


48 


10, 


.591 


+33 


43 


27. 


,13 


0, 


.356 


0, 


,64 


1, 


,18 


140 + 11.7 


2, 


,03 


44, 


,4 


± 


2, 


.3 


12, 


.64 


8.47 


2, 


.75 


FF 0758+2851 


07 


58 


45, 


.956 


+28 


51 


32. 


,76 


0, 


.126 


0, 


.62 


0, 


.92 


87.5 ± 9.52 


3, 


,72 


40, 


.0 


± 


2, 


.2 


11 


.58 


7.64 


2, 


.42 


FF 0804+3919 


08 


04 


07, 


.399 


+39 


19 


27. 


,63 


0, 


.164 


0. 
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0. 
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5. 
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,2 


± 


2, 


.7 
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1, 
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.0 
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8.28 


2, 


.70 


FF 0826+3042 


08 


26 


11, 


.644 


+30 


42 
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35, 


,7 
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1, 


.6 
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8.46 
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0. 


,51 
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0. 
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50 
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.220 


0, 


,24 


0, 
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8, 


,51 
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,1 


± 


1, 


.5 
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.93 


8.28 


1, 


.88 
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07 


42, 


.264 


+39 


31 


49. 


,47 


0, 


.224 


0, 
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0, 


,47 
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1, 


,24 


41, 


,6 


± 


2, 


.2 
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2, 


.56 


FF 1016+3951 
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16 
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.616 
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51 
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0, 


.307 


0, 
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0, 


,39 
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2, 


,10 


43, 


,3 


± 


2, 


.2 


12 


.01 
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.25 
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18 


34, 


.539 


+36 


49 
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0, 


.490 


0. 
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0. 


,56 
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9. 


,18 
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,2 


± 


1, 


.9 


12 
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8.42 


1, 


.72 
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42 
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31 
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0. 
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0, 
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6, 
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,8 
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.5 
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06 
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.716 
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01 


46. 


,39 
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0, 


,19 


0, 


,64 


62.5 + 16.2 


13.7 


57, 


,9 


± 


2, 


.5 
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.09 


8.40 


1, 


.44 


FF 1242+2905 


12 


42 


32, 


.497 


+29 


05 


14. 


,75 


0, 


.260 


0. 


,19 


0. 


,50 


85.3 + 26.5 


1. 


,34 


36. 


,6 


± 


1, 


.6 


11, 


.95 


8.21 


2, 


.55 


FF 1412+3014 


14 


12 


24, 


.952 


+30 


14 


09. 


,78 


0, 


.257 


0, 


,15 


0, 


,53 


63.3 ± 14.2 


2, 


,06 


36, 


,8 


± 


1, 


.5 
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.86 


8.11 


2, 


.29 


FF 1456+3337 
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56 
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.427 
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37 


09. 
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0, 


.443 


0. 
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0. 
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1. 
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,4 
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.1 
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14 
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29 


42. 
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0, 
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(22.5 + 10.6) 
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.052 
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42 
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0, 


.926 
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,26 
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52.3 + 13.4 
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,89 


62, 


,6 
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2, 


.8 


13 


.12 
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.81 
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14 
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.105 
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34 
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2, 
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17 


07 
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25 
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0, 
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0, 


,13 
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2, 
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.4 
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.08 


8.44 


2, 


.31 


FF 1713+3843 


17 


13 


46, 


.085 


+38 


43 


04. 


,77 


0, 


.171 


0, 


,37 


0, 


,79 


72.1 + 21.2 


1, 


,66 


39, 


,1 


± 


1, 


.9 


11 


.72 


7.82 


2, 


.62 


FF 2131-0141 


21 


31 


53, 


.490 


-01 


41 


43. 


,35 


0, 


.730 


0. 


,11 


0. 


,47 


64.9 + 21.2 


2, 


,79 


48, 


,1 


± 


1, 


.9 


12, 


.82 


8.48 


2, 


.07 


FF 2136-0112 


21 


36 


34, 


.229 


-01 


12 


08. 


,38 


0, 


.210 


0. 


,29 


1. 


,11 


69.5 + 25.8 


3. 


,20 


38, 


,9 


± 


1, 


.7 


11 


.86 


7.98 


2, 


.29 


FF 2200+0108 


22 


00 


51, 


.859 


+01 


08 


27. 


,08 


0, 


.164 


0. 


,20 


0. 


,72 


( ••• ±41.9) 


2. 


,82 






















FF 2216+0058 


22 


16 


02, 


.721 


+00 


58 


10. 


,65 


0, 


.212 


0. 


,51 


0. 


,85 


48.0 ± 16.4 


1, 


.31 


46. 


,8 


± 


2, 


.7 


11 


.96 


7.68 


2, 


.76 



Table 1 — Continued 



Source Name 


RA 


Dec 


z 




SlOO/Lxm^ 




Sl.4GHz 


Td 


log Lfir 


log Md 


5 




J2000 


J2000 




Jy 


Jy 


mJy 


mJy 


K 


L© 


M© 




FF 2221-0042 


22 21 26.066 


-00 42 39.08 


0.189 


0.20 


0.37 


42.8 ± 12.8 


0.66 


39.1 ±2.0 


11.54 


7.65 


2.75 


FF 2330-0025 


23 30 34.920 


-00 25 03.98 


0.252 


0.25 


0.72 


137 ±33.8 


1.65 


35.2 ±1.5 


12.08 


8.42 


2.62 


FF 2352-0015 


23 52 53.171 


-00 15 24.69 


0.227 


0.19 


0.61 


(88.8 ± 48.6) 


0.59 











^The template amplitude of the noise- weighted mean scan (1003), given by IRAS Scan Processing and Integration (SCANPI), Version 5.0. 

^The template amplitude of the noise-weighted mean scan (1003), given by SCANPI, Version 5.0, except for FF 0123+0114 and FF 2200+0108, for 
which the 100 jim fluxes are not available from SCANPI and thus values in Stanford et al. (2000) are adopted. 



Note. — Three additional far-IR/submm flux measurements are available from the literature - FF 1532+3242: Sgoi_^rn= 478 mJy, Si80f^m= 397 mJy 
(Verma et al. 2002); FF 1614+3234: S850fxm= 8.47 mJy (Farrah et al. 2002). 

Note. — /3 =1.5 is assumed in the SED fltting. 
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Fig. 1. — The redshift distribution of the 36 FF sources selected for SHARC-II observations. 
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Fig. 2. — 350 /xm sigma-to-noise maps of the 28 FF sources detected by SHARC-II. 
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Fig. 2. — Continued 
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Fig. 3. — 350 /xm sigma-to-noise maps of the eight FF sources not detected by SHARC-II. 
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Fig. 4. — Far-IR/submm SED fits of the SHARC-II detected intermediate-redshift ULIRGs. 
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4. — Continued 



- 22 - 




Td / K 

Fig. 5. — Histogram of Ta values derived for our intermediate-redshift ULIRG sample. 



-23- 




Log(Lp,R/y 



Fig. 6. — a): The far-IR/radio correlation for our intermediate-redshift ULIRG sample. The 
five sources which have been robustly identified as either AGN or starburst (SB) dominated 
systems have been highlighted. The solid line represents a linear fit to the 21 sources with 
q > 2.0, and the dotted lines mark the ±0.2 scatter in q. The dashed line is the local far- 
IR/radio correlation given by Condon et al. (1991). b): The distribution of q values as a 
function of the far-IR luminosity. The solid line shows the sample median for all 28 sources, 
and the dotted lines mark the ±0.4 scatter. The dashed line shows the median of the 21 
sources with g > 2.0. 
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Fig. 7. — LpiR and derived for luminous dusty galaxy samples at distinct redshifts. 
The dashed lines correspond to Lpm-Td relations given by Eq. f fTOl ). assuming f] =1.5 and 
= 1.875 (kg/m^)~-^, for dust mass spanning over three orders of magnitude = 
10^"^° Mq (from top to bottom). 



